Topological Schrodinger cats: Non-local quantum superpositions of topological defects 
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Topological defects (such as monopoles, vortex lines, or domain walls) mark locations where 
disparate cho ices of a broken symmetry vacuum elsewhere in the system lead to irreconcilable 
differences^! They are energetically costly (the energy density in their core reaches that of the 
prior symmetric vacuum) but topologically stable (the whole manifold would have to be rearranged 
to get rid of the defect). We show how, in a paradigmatic model of a quantum phase transition, 
a topological defect can be put in a non-local superposition, so that - in a region large compared 
to the size of its core - the order parameter of the system is "undecided" by being in a quantum 
superposition of conflicting choices of the broken symmetry. We demonstrate how to exhibit such a 
"Schrodinger kink" by devising a version of a double-slit experiment suitable for topological defects. 
Coherence detectable in such experiments will be suppressed as a consequence of interaction with 
the environment. We analyze environment-induced decoherence and discuss its role in symmetry 
breaking. 



Topological defects are the epitome of locality. An ex- 
ample of a defect occurs in the quantum Ising model 
where a lattice of spins interact ferromagnetically with 
their nearest neighbors, i.e., the Hamiltonian contains 
an interaction of the form — c^ct^ +1 for the n th spin on 
the lattice. The entire collection of spins achieves its low- 
est energy when they are all aligned. However, there are 
two choices for this lowest energy state, | • • • tttt ' ' ' ) or 
| ■ ■ • J444- ' ' ' ) ■ Both are energetically identical but each 
of them breaks the symmetry of the Hamiltonian, which 
has no preference between "up" and "down" . 

A configuration that includes topological defects could 
arise, for example, in a driven phase transition, i.e., "a 
quench" . In the Ising model, a quench can be induced 
when the transverse field strength g is changed in the 
Hamiltonian 

n 

When the external field is strong enough (g ^> 1), it 
"wins" and all spins align with a x . A decrease in g, 
however, will lead to a phase transition at g = 1 with the 
symmetry breaking term —cr^a^ +1 trying to align all the 
spins in their z-direction. 

The choice of whether the spins will point up or down 
is made locally. As a result of causality and the finite 
speed at which signals propagate, the size of these do- 
mains will be determined by rate of change of g^^. This 
will lead to configurations, such as | • ■ ■ tt44- where 
the topological defect marks the location of a "switch" 
from one broken symmetry ground state to the other. 
In this one-dimensional example, the defect is a kink (a 
domain wall), but, in a sense, all the topological defects 
(monopoles, vortex lines, etc.) that exist in higher di- 
mensions "look the same". A discussion of the conse- 
quences of this process for the post-transition state is 
beyond the scope of this work (but see, e.g., Refs. ISr ITTj) . 



We only note that the density of the resulting kinks will 
be finite. 

The quantum Ising model quenches a state accord- 
ing to Hamiltonian ([!]) so, when g is changed in time, 
superpositions of different locations of a kink (e.g., 

a\ ■ ■ ■ mm ■ ■ ■ ■ ■ ■ mm ■ • • • • • ttttu • • • » 

are allowed, and, indeed, inevitabl e! 10 * 11 !. Spreading of a 
localized kink will come about as a consequence of the 
kinetic term, gcr^, in equation 0. A convenient setting 
to investigate its effect is offered by a slight modifica- 
tion of the original Hamiltonian, so that the Ising in- 
teraction between a selected pair of sites, no, no + 1, is 
— (1 — w) o-f l0 <7„ 0+1 with w > 0, which differs somewhat 
from the uniform coupling of —a^a^ +1 . This difference 
means that the kink is energetically less expensive when 
localized between those two selected sites. The decrease 
in the coupling constant by w creates a local "potential 
well" that binds this kink. On the other hand, the ki- 
netic term will delocalize the kink so that the quantum 
wavefunction of the kink will have the form 

^oce-^l, (2) 

where i\) n is the amplitude for the kink to be on the link 
n between sites n and n + 1, no is the location of the po- 
tential well, and 70 = sinh -1 (w/g) is the inverse decay 
length of the wavepacket. One can imagine measure- 
ments that will reveal such a non-local wavepacket - a 
kink in a superposition of many locations. We empha- 
size that the half-width of this wavepacket is not the size 
of the kink: The kink is a local object with a size given 
by the healing length that - in the quantum Ising model 
far away from the critical point - is given by the lattice 
spacing of the neighboring spins. The spread in equation 
([2]) represents a superposition of many possible locations 
of the kink, which is bound to the weak link between sites 
no and no + 1, but nevertheless has some spatial extent. 
A tell-tale signature of quantum coherence is an inter- 
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Figure 1: A Schrodinger kink in a quantum Ising chain, (a) A topological defect in a non-local superposition and (b) The 
analogue of a double-slit experiment. A double-well potential (left) is used to create a topological defect, such as a domain wall, 
superimposed in two locations (here, defect is represented by its probability distribution in space) . To carry out the double slit 
experiment, the two potential wells are eliminated, allowing the defect to move. In isolation, the two wavepackets emerging 
from the "slits" interfere creating fringes of high and low probability for the location of the kink. However, when the system 
interacts with the environment, such a superposition will decohere at rate proportional to the distance L - the unzipped part 
of the chain shown in (a), which corresponds to the size of the Schrodinger kink - resulting in a classical distribution for the 
defect. 



ference pattern. To see whether a defect can interfere 
with itself, we devise an analogue of the double-slit ex- 
periment (see Figure [T] for a schematic). To this end, the 
local value of the Ising coupling can be depressed by w 
at two locations separated by L sites. The two links that 
bind the same kink initially are analogues of the slits in 
the double-slit experiment. To achieve a situation where 
the kink is "both here and there" one can start it on one 
of the two binding sites and evolve so that tunneling of 
the kink will result in the state 



ttt 11 



(3) 



where we have ignored the spread of the wavepacket for 
simplicity. This is also illustrated in Figure [lja), where 
the superposition of the kink forms a "scar" in the ori- 
entation of the spins, "unzipping" the L-sized region of 
the spin chain in Fig. 1. 

This is not the only way to create a "Schrodinger kink" 
state. As we shall see later, one can also start with a kink 
localized state on a single "weak link". If the kink is 
released symmetrically, it will travel both left and right. 

The double-slit experiment can be now conducted 
starting from this initial configuration by eliminating the 
two weak links. The kink is then no longer bound by 
the "double-well" potential: The Schrodinger kink prop- 
agates in accord with the Schrodinger equation and the 
two components of the wavepacket run into each other. 
If coherence was properly preserved between them, this 
will lead to an interference pattern. 



Preserving phase coherence is crucial if the resulting 
interference pattern is to be seen in experiments. Triv- 
ial reasons for loss of coherence - such as an imprecise 
implementation of the two weak links - will have to be 
eliminated. However, the fundamental reason for the loss 
of coherence is environment- induced decoherence^HUl. It 
is important to understand its causes and its nature, as it 
is not just an impediment to creating superpositions de- 
scribed above, but the prevailing reason why topological 
defects we encounter are always localized. 

Decoherence is caused by the interaction of individual 
spins with the environment £. For instance, an individ- 
ual spin | f) can leave a different imprint on £ than the 
spin | I), i.e., (a| t)+/3| i))\£o) a| t)|£ t ) + P\ IM)- 
The overlap {£^\£^) determines the decoherence factor, 
with (£^\£i) = corresponding to the complete loss of 
coherence. The decoherence factor controls size of the 
off-diagonal terms in the density matrix. When they 
disappear, quantum coherence between | f) and | I) is 

Returning to our Schrodinger kink, we note that when 
two locations are separated by L spins, equation the 
decoherence process will take place simultaneously in all 
L spins. Assuming that each spin leaves its own imprint 
will lead to a decoherence factor that scales as (£^-\£i) L , 
where L is the number of unzipped spins - the spatial ex- 
tent of the superposition of the "Schrodinger kink" . This 
exponential scaling with the extent of the superposition 
is intuitively obvious: We do not have to assume any 
specific model for the decoherence. All that is needed 
is the familiar assumption (see, e.g., modeling of errors 
in quantum error correctiouISl) that individual spins (or 
local regions) affect the environment individually. 
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Figure 2: Interference patterns after a Schrodinger kink is re- 
leased, (a) Time evolution of the Schrodinger kink for L = 50 
and 2w = 0.3g results in an interference pattern (highlighted 
in red at the final time), (b) Interference pattern in the long 
time limit (gt = 1000) for L = 100 and 2w = 0.3g. The exact 
data is shown as black crosses and the expression, equation 
Q, is shown as the blue line, (c) Interference pattern in 
the long time limit (gt — 1000) for L — 50 and 2w — 0.3g. 
The exact data is shown as black crosses and equation Q is 
shown as the blue line. As the initial superimposed kink is 
shortened, the interference pattern diminishes. 
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Figure 3: A Schrodinger kink evolved in the presence of deco- 
herence. (a) Time evolution for strong decoherence (F = g/2), 
L — 50, and 2w = 0.3g, where the Schrodinger kink evolves 
into a (b) Gaussian mixture of locations (highlighted in red 
in (a)) at a later time, gt — 100 (the black crosses are the 
exact data and the blue curve is the solution to the diffusion 
equation), (c) Under weak decoherence (r = gl0~ 3 /8), the 
superposition is still visible at intermediate times, but the 
decoherence smooths out the fringes (the black crosses are 
the exact data and the blue curve is the pure state solution 
convoluted with a Lorentzian). 



strated explicitly by the form of the fringes, 



This assumption suffices to show that the decoherence 
rate is proportional to the "size", L, of the superposition. 
That is, decoherence time is T& ec ~ 1/L. This conclusion 
is confirmed by calculations employing a master equation 
(see the Supplemental Information). One can generalize 
this intuition to superpositions of topological defects in 
higher dimensions by noting that it is the volume of the 
system that is the superposition, i.e., that is suspended 
in indecision between two broken symmetry vacua, that 
is responsible for the decoherence rate. 

We now return to the Ising model. In the absence of 
decoherence, a well-defined interference pattern develops 
when the two initial components of the wavepacket run 
into each other (Fig. [5]). As the distance, L, between 
these two components is increased, more fringes become 
visible (i.e., decreasing L decreases the magnitude of the 
outer fringes), as seen by Figs. [2^b,c). This is demon- 



Pn (t) = \lpn (*)l °C 



1 + COS 



(n-n -L/2)L 



1 



(n— n — L/2) 2 
(2709*? 



(4) 



This interference pattern is analogous to the double-slit 
experiment. In the double-slit experiment the distance 
between fringes is \D / L, where L is the distance between 
the slits, A is the wavelength of light, and D is the dis- 
tance to the screen. For kinks, the distance traveled is 
D = Vkt ~ 2gkt — 2g (2tt / X) t, where Vk — duj/dk ss 2gk 
is the group velocity of the kink for slow kinks that are 
generated when w -C 1, as we assume in deriving Eq. Q 
(the figures use Vk = 2gsink, demonstrating that this 
assumption is sound). This linear approximation yields 
the fringe distance for kink interference XD/L = Awgt/L, 
as given in equation Q. The height of the second peak 

relative to the first is 1/ [l + 47t 2 /L 2 7q] 2 . In effect, the 
linear approximation of Vk = 2g sin k entitles one to think 
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Figure 4: A single kink evolved on a finite lattice with and without decoherence. (Left) Time evolution without decoherence 
on an L = 201 site lattice and 2w = 0.5g. The kink travels outward on the lattice and after reflecting off the boundaries, 
it starts to interfere with itself. The leftmost panels show the self interference at two particular times. (Right) The same 
simulation except in the presence of weak decoherence (F = 3g x 10 -6 ). The kink can still interfere with itself, but eventually 
the decoherence will set in. In the supplemental information, we show movies of the development of this interference pattern 
for different values of w. 



of kinks in a double slit experiment as one would think 
of, say, electrons or neutrons. 

Decoherence suppresses interference. Figure p[a,b) 
shows the evolution of the kink probability density (that 
determines the probability of finding a kink at a certain 
location) under strong decoherence. In presence of deco- 
herence it evolves into a Gaussian form seen in Fig. [3jb) . 
Even weak decoherence will attenuate interference. The 
decoherence time depends on L: Coherence disappears 
exponentially fast in L, the separation between compo- 
nents of the wavepacket. The evolution of the superpo- 
sition under weak decoherence is shown in Fig. (3^c). 

The double slit experiment for kinks we have just dis- 
cussed has the advantage of a straightforward interpre- 
tation. The use of the coherent bi-local Schrodinger kink 
that is relatively weakly bound to the two sites allows 
one to linearize the dispersion relation and use = 2gk 
for the kink velocity. This leads in turn to the simple 
form of the interference pattern, equation Q. 

However, this ease of interpretation may come at the 
price of difficult implementation. In particular, prepar- 
ing the initial bi-local kink wavepacket and maintaining 
coherence between its two pieces of the will be a chal- 
lenge. A different version of Schrodinger kink that should 
be easier to implement is therefore illustrated in Fig. [4] 
Here the kink is initially bound to a specific link along the 
Ising chain. This kink trap is instantaneously turned off, 
which results in a coherent spreading of the wavepacket 



with a superposition of velocities and in both directions 
on the Ising chain. 

As before, we are not satisfied with just creating a 
Schrodinger kink. One should confirm that quantum co- 
herence is present. In Fig. [4] this is accomplished by 
reflecting the spreading wavepacket from the ends of the 
Ising chain. The time-evolving interference pattern is 
now more complicated than before, but - in absence 
of decoherence - it clearly exhibits quantum coherence 
that can be predicted by suitably "folding" the kink's 
wavefunction upon itself. Decoherence (as expected) sup- 
presses interference fringes over time: As with double slit 
analog, the decoherence strength will have to decrease as 
~ 1/L, where L is the support of the wavepacket, if co- 
herence - and, hence, interference - is to survive. 

The interference pattern in Fig. 4 bears the imprint of 
the dispersion relation on a lattice, u> = ~2g cos k: When 
released, a tightly bound kink turns into a wavepacket 
that propagates as a Bessel function, p n (t) = \J n (2gt)\ 2 . 
Our kinks of both Fig. 4 and especially of Fig. 2 are 
relatively weakly bound. Therefore, the Bessel oscillation 
is suppressed - smoothed out by the non-local nature of 
the wavepacket that eliminates large k contributions (see 
Supplementary Information). Nevertheless, small-scale 
jaggedness of the interference pattern visible in Fig. 4 
(where w is larger than in Fig. 2, and, therefore, the 
kink starts more tightly localized) is a remnant of these 
Bessel oscillations. 
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An obvious application of this observation is the "col- 
lapse" of the superposition of the broken symmetry vacua 
after a phase transition. For example, in the case of the 
quantum Ising model, the ferromagnetic ground state is 
a superposition of | t • • • f) and | I • • • l) . The total 
number of spins, N, in the Ising chain is the size of the 
superposition (e.g., | t " * ' t) + I 4- " ' but this sym- 
metric superposition will become a mixture of the two 
obvious broken symmetry states in a very short time, 
Tdec ~ 1/N. This is a simple and compelling explanation 
of the symmetry breaking that occurs whenever a phase 
transition starting from a symmetric vacuum takes place. 

The effectiveness of decoherence in localizing topolog- 
ical defects provides novel insights into symmetry break- 
ing dynamics. A phase transition leads to a single quasi- 
classical configuration pockmarked with topological de- 
fects: Decoherence is the final step in this process. It 
leads to the "collapse of the wavepacket" that initially 
contains all the possible broken-symmetry configurations. 

In the end, only a simple quasi-classical configura- 
tion is in evidence. This is essentially the same course 
of events that takes place in quantum measurement!^, 



where breaking of the unitary symmetry allows for a "col- 
lapse of the wavepacket" . Given the recent progress in 
emulating quantum Ising and other models^"^!], we be- 
lieve that experimental tests of "Schrddinger kinks" are 
within reach. Moreover, although models of supcrflu- 
ids do not allow one to develop and analyze microscopic 
quantum theories of double-slit experiments for the rele- 
vant topological defects with the detail we have presented 
above for kinks, experiments involving the generation 
and manipulation of such defects in, for example, gaseous 
Bose-Einstein condensates allow one to hope that testing 
of 'topological Schrddinger cats' involving vortex lines or 
solitons may also be possible. 
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Supplemental Information 
Interference fringes 

We are interested in the one-kink subspace of the 
Hamiltonian (1) spanned by states 



\n) = | tt ••• tnWl ••• W>, 



(5) 



which have a kink on the link between sites n, n + 1. We 



can confine to this subspace when j < 1 



and the 



magnetic field g is too weak to mix with the subspaces of 
3, 5, .. kinks. When there is only one weaker link, then a 
stationary state \tp) = X^VViln) satisfies a Schrodinger 
equation 



Elp„ = -g (tpn+l + lpn-l) - 2wS no ^ n . 



(6) 



The magnetic field g provides a hopping term between 
nearest neighbor links and the weaker link of strength 1 — 
w results in a trapping potential of strength 2w localized 
on the link n . This potential has one localized bound 
state 

ip n = cosh (70) e~ 7 °l n ~" D l, (7) 

where 70 = sinh -1 (w / g) , with energy E = 
— 2gcosh(7o) = —2\Jg 2 + w 2 . 

When there are two weaker links, no and no + L, the 
Schrodinger equation 

Elpn = -g(lp n +l + 1pn-l)- 2w(S nQ . n + 8 no+L ,n) VVs (8) 

has two bound states 

oc e" 71 ™-™ 01 T e'^ n - na - L ^ (9) 



with energy E = — 2g cosh(7) and 7 being one of the two 
solutions of the equation 



g sinh (7) 
w 



(10) 



The greater of the two 7's corresponds to the symmetric 
ground state ip + . 

We are interested in the "tight-binding" regime where 
the wave packets in equation (J9J) overlap weakly, i.e., 
e -7_L ^- 1 ln this regime 7 70 in both bound states 
?/' =F , their energies are E ss Eq ± w, respectively, where 
the gap 



2uj 



Aw 1 



-2 l0 L 



Vi 2 



(ii) 



is relatively small, 2oj <C Eq. In this regime, we can 
initially prepare the ground state of a single well, equa- 
tion which, in the basis reads ip + + tp~ ■ Af- 
ter this preparation, we can either switch on the second 
well suddenly or turn it on adiabatically. An alternative 



preparation is to split a single well adiabatically and sym- 
metrically into two wells. With a real-time preparation, 
the initial state evolves into ^e" 1 " 2 "' +ip~e~' tult and the 
state becomes an equal superposition of the left and right 
wells, e" 701 ™"™ 01 + i e -io\n-n -L\ ) at the ear ii est tmie 

ujt = jTT. In the adiabatic case, the preparation ends 
in the ground state ip + of the double well. Since both 
cases require roughly the same time ~ we opt for 
the more robust adiabatic preparation of ip + . 

Once the state ip + has been prepared, we switch off the 
double-well potential at t = 0, w — > 0, to let the wave 
function freely disperse with just the hopping term, 



d 1 



-9 OA 



n+l 



(12) 



When 2(/7q t ^> 1, the probability distribution develops 
an interference pattern 



p n (t) = \lp n (t)| CX 



1 



(w-n -L/2)£ 
2gt 



1 



(rt-Ko-L/2) 2 
(270S*) 2 



(13) 



Here the distance between fringes is Airgt/L and the 
width of the Lorentzian-squared envelope is 2jogt. In 
the tight-binding regime, where 70-^ 3> 1, we obtain a 
large number of fringes. Equation ( 13 1 is the result we 



plot with the exact solution within the main text. 

If, on the other hand, one starts from a single kink 
localized on link 0, ip n (t = 0) = 5 7lt Q, it will evolve into 



dke 2l9t 



cos k ^ikn 



J n (2gt), (14) 



where J n is a Bessel function. This will also give inter- 
ference, but not of the form in equation (13). We are 



interested in the regime where these "intrinsic" oscilla- 
tions due to the lattice are negligible. That is, if we start 
with a kink wavepacket with some spread, 70 from above, 
then 



dke 2i3tcosk e ikn f[k/j ], (15) 



where / will give the Fourier transform of the initial 
wavepacket that cuts off large k compared to 70. When 
the kink is initially tightly bound tightly (large w and 
70), then / is a constant for — tt < k < tt and we ob- 
tain ~ J n (2gt). When w is weak and 70 <S 1, then 
/ cuts off |fc| greater than 70 and we can approximate 
cos A; rj 1 — k 2 /2, i.e., we can linearize the group velocity 
and obtain equation (13). The three movies show the 



development of the interference pattern (of Fig. 4) for 
different values of w. 



Decoherence 

Once the double-well potential has been switched off 
there are no energy gaps 2uj or Eq to protect against even 
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weak decoherence. Under local, Markovian dephasing, 
the state, described by a density matrix pit), evolves 
according to a master equation 



dt 



p(t) = -i[H ,p(t)] 



4^ 



K,K, P (t)]], (i6) 



where Hq is the hopping Hamiltonian that gives rise 
to the evolution in equation (12 1 and the initial state 
is p(0) = \ip + )(4> + \- In the position representation, 
P = Em n Pm,n\'m){n\, the master equation reads 



dpm,r 

dt 



-W {Pm+l,n + Pm—l,n ~ Pm,n+1 ~ Pm,n-l) 

- T\m- n\ p m ,n- (17) 



We can consider regimes of strong or weak decoherence 
when either r > j or T < g respectively. 



Strong decoherence 

When we temporarily set 3 = 0, the off-diagonal ma- 
trix elements decay like p m , n (t) — p(0) exp(— F\m — n\t). 
Therefore, we can assume that deep in this regime only 
the diagonal elements, p n ^ n = p n , and near-diagonal ele- 
ments, p n .n+i = z n and p n +i,n = z*, are non-zero. All 
other elements, even if they are non-zero initially, quickly 
become negligible. The master equation then simplifies 
to the set of equations 



dp n 
dt 



g (Sn-l - S n ) 



and 



ds n 
dt 



-Ts n - 2g (pn+l ~ Pn) , 



(18) 



(19) 



where s n — —i (z n — z*). Since V ^> g, the slave field s n 
can be adiabatically eliminated and we obtain a lattice 
diffusion equation 



dt 

with the diffusion constant 



D 



._, = D (p n+1 - 2p n + Pn-t) 



2g 2 



(20) 



(21) 



When 2Dt >• 1, the initial distribution p n = \ipt,\ 2 
spreads into 



pp >9) (t) = Ei^ 



i 



-e iDt 



(22) 



without any interference fringes. 

Weak decoherence 



In case of weak decoherence, r <C g, the initially 
smooth wavefunction does not get localized in space so 
we can make a long wavelength approximation (LWA) 
where the lattice site numbers m, n are continuous coor- 
dinates. In quasimomentum representation, 



Pp,q( f ) = X! Pm < n ex P[*(P TO _ 1 n ) 



(23) 



after going to the interaction picture and making the 
LWA, 



P P ,q{t) = p p ,q(t) exp[2igt(cosp- cos q)} 

« p p , q {t)exp[-igt(p 2 -q 2 )}, (24) 



the master equation (17) becomes 



dp Pl q 
dt 



dk 
2~i 



e - 2tg t( P - q )k ~ p+k q+k / dm \ m \ e 

J — OO 

(25) 

Its exact solution gives a probability distribution in space 



-ikm 



Pn r<<9 Ht) = 1 



\ ^ 



i(t) 



it ^ l 2 (t)+m 2 



Pn-\-m 

(t), (26) 



where l{t) = gTt 2 and p m (t) is a probability distribution 
in the absence of decoherence, i.e., the interference fringes 
in equation (4). 



The Lorentzian convolution |26|) coarse-grains the 
fringes in p m (t) on the scale l(t). l(t) becomes greater 
than the distance between fringes at ^dec — 47r/rX. This 
is the (not quite unexpected) decoherence time when the 
environment can distinguish between the f and -J, mag- 
netization of the L sites between the potential wells. 



The convolution (26 1 tends to a Lorentzian when l(t) is 
much greater than the width of the envelope in equation 
(4) or, equivalently, Tt 3> ^70 ■ The Lorentzian is differ- 
ent (wider) than the Lorentzian-squared envelope in the 
decoherence-free fringes (4) and the asymptotic Gaussian 
(22) in the strong decoherence limit. 



